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Abstract 

The phase diagram of the O(n) model, in particular the special case 
n = 0, is studied by means of transfer- matrix calculations on the loop 
representation of the O(n) model. The model is defined on the square 
lattice; the loops are allowed to collide at the lattice vertices, but not 
to intersect. The loop model contains three variable parameters that 
determine the loop density or temperature, the energy of a bend in a 
loop, and the interaction energy of colliding loop segments. A finite- 
size analysis of the transfer-matrix results yields the phase diagram 
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in a special plane of the parameter space. These results confirm the 
existence of a multicritical point and an Ising-like critical line in the 
low-temperature O(n) phase. 



Keywords: O(n) model; Polymers; Phase diagram. 

2 



1 Introduction 



The 0(n) model is originally defined as a model of n-component spins, with 
O(n) symmetric interactions. This model appears to be meaningful only 
when n is a positive integer. However, a mapping [p]] has been found on a 
model of closed loops in which n appears as the Boltzmann weight associated 
with each loop: hence n has become a continuously variable parameter. The 
critical point, which is characterized by e.g. the divergence of the susceptibil- 
ity in the spin representation, corresponds with a diverging loop size in the 
loop model. The latter model, in particular in two dimensions, is accessible 
by exact and numerical analyses, and is thus also relevant for the exploration 
of the O(n) spin model. Indeed, a wealth of results have been found for the 
O(n) loop model. For instance, an exact critical line was found for an O(n) 
model on the honeycomb lattice, as well as the associated thermal and mag- 
netic exponents p], |3], |], ||. The critical line covers the interval —2 < n < 2 
of the loop model and thus applies to n = 1 and n = 2 in the spin model. 
Another special case of the O(n) loop model is n = 0, where the weight of a 
loop vanishes. This does not mean that the system is empty: the O(0) model 
may be formulated in terms of the derivative of the O(n) partition function 
to n, taken at n = so that one loop remains, finite or infinite. Thus the non- 
intersecting O(0) loop model can be applied to describe self-avoiding walks 
(SAW's) or polymers. Furthermore the n = —2 case corresponds to random 
walks without excluded- volume interactions || [7|. 

More complicated phase behavior is observed when attractive interactions 
between the loop segments are introduced. As long as these attractions are 
sufficiently small, the O(0) ordering transition remains in the same univer- 
sality class. But when the interactions reach a certain threshold a tricritical 
point appears, called the theta point [§], and at even stronger attractions 
the phase transition turns first order. While exact information is scarce in 
the temperature versus attraction phase diagram, a numerical analysis || 
confirms the picture sketched here. 

This tricritical theta point was first found || for a loop model with va- 
cancies on the honeycomb lattice. The same type of transition, describing a 
collapsing polymer, was found for the square lattice O(0) loop model, namely 
on branch of Ref . JTOj] . The latter model did, in addition to a temperature 



variable and an attraction between loop segments, also contain a parameter 
that distinguishes between straight loop segments and 90 degree bends. The 
results obtained for the square lattice included a generalization |Ttl of the 
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theta transition to values Remarkably, the enlarged parameter space 

of three parameters (not counting n) appears to contain also a different type 
of collapse transition p0| , |TT|j for n = 0. Like the theta point found on the 
square lattice, it generalizes to other values of n (branch 3 of Ref. fllPp . 

To understand this second collapse mechanism, it may help to identify 
straight loop segments in terms of Ising-like degrees of freedom in the O(n) 
model |l(J |9f . When the straight loop segments are suppressed, Ising-like or- 
der appears. Thus the newly observed multicritical point has been identified 
as the point along the O(n) critical line where the Ising degrees of freedom 
become critical. These new effects are a consequence of the underlying lattice 
structure and are thus not generally expected in a model of free polymers. 
However, they could be present in the case of a polymer influenced by a 
periodic potential caused by a crystalline substrate. 

The freezing of these Ising degrees of freedom does not only occur on the 
line of phase transitions, but is also expected in the low-temperature 0(n) 
phase. Thus an Ising-like transition would extend from the multicritical 
point into the dense phase of the loop model. Indeed, the results reported in 
Ref. [plj include a branch of critical points where the critical exponents reflect 



both the Ising and the low-temperature O(n) universal critical behavior. 

These observations have been interpreted in terms of a conjectured and 
qualitative phase diagram [III. In the present paper we confirm the validity 



of this phase diagram, including the predicted Ising line, by means of finite- 
size scaling and a transfer-matrix technique. Our investigations concern the 



n = case of the O(n) model on the square lattice as defined in Ref. JlO 
We briefly recall the definition of the partition function 

Z loop = £nVV^ (1) 
c 

which depends on three vertex weights x, y, and z, associated with a 90- 
degree bend of a loop segment, a straight loop segment, and a vertex visited 
by two 90-degree bends: a collision. Intersections are forbidden and a weight 
1 is assigned to an empty vertex (not visited by a loop). The exponents N x , 
N y and N z denote the total number of vertices of the type indicated. The 
number of closed loops is iVj, and thus, after differentiation with respect to n, 
only loop configurations C with jVj = 1 contribute to the partition function 
of the n = model. When the weight y of the straight loop segments is 
set to 0, we freeze the Ising-like degrees of freedom mentioned above. Then 
the model reduces to that investigated in Ref. |§ and displays a theta point. 
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The special character of the y = case also becomes apparent from the fact 
that the model, at least for even system sizes, is symmetric with respect to a 
change of sign of y. The point y = is special on the basis of this symmetry. 
In this paper we explore nonzero values of y and restrict ourselves to the 
x = y plane. The weight x — y plays the role of a temperature-like parameter 
(at least when z is small) and drives the SAW-type transition between the 
vacuum and the dense phase. The role of the parameter z is to adjust the 
attractive interaction between loop segments. 

The rest of this paper is composed as follows. In Section 2 we summarize 
the numerical procedures and in Section 3 we present the results and discuss 
the phase diagram thus obtained. 



2 Numerical procedures 

We consider the square O(0) model, wrapped on an infinitely long cylinder of 
circumference L, the finite size of the system, which is measured in elemen- 
tary lattice spacings. One can derive quantities of interest, such as the free 
energy and correlation lengths, from the eigenvalues of the transfer matrix. 
For this purpose we make use of the construction of the 0(n) loop model 
transfer matrix given in Ref. [ID"|. 

The relevant eigenvalues of the transfer matrix are obtained by a method 
explained e.g. in Ref. |T2|. This method is based on iterative multiplication 



of a vector by the transfer matrix. Even when n = this procedure gener- 
ates non-empty configurations containing a part of a loop. Closure of such 
a loop leads to a zero contribution to the resulting vector, but there will 
always be contributions due to configurations in which the loop is not (yet) 
closed. These configurations contribute to the eigenvalues of the transfer 
matrix. Therefore the transfer matrix allows the exploration of non-empty 
configurations of the O(0) model. 

The 0(n) spin-spin correlation function can simply be translated into 
the language of the 0(n) loop model, as the ratio Z'/Z, where Z is the 
partition sum of the loop model and Z' a similar sum, but restricted to 
loop configurations that include a non-closed loop segment connecting the 
two correlated spins. For spins separated along the length direction of the 
cylinder, we may obtain Z' from the 'odd transfer matrix' that describes 



configurations which include such an additional single segment [TO]. The 



partition zum Z is obtained from the 'even transfer matrix' that allows only 
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configurations in which loop segments are pairwise connected at the end of 
the cylinder. This is a convenient method to derive the magnetic correlation 



length and its associated scaling dimension. For further details see Ref. fLO 

The relation between the magnetic correlation length the largest eigen- 
value of the even transfer matrix and the largest eigenvalue A^ of the 
odd transfer matrix is 

aW = l/ln(Ai 0) /Ai 1) ) (2) 

The finite-size scaling behavior of £/, and the principle of the 'phenomeno- 



logical renormalization' [13] of ^ provide a method to explore the phase dia- 
gram of the model. For this purpose we define the 'scaled gap' as X h (t, L) = 
L / 2n £h{t, L). Apart from the finite size L, it depends on a temperature-like 
scaling field t. More than one scaling field can play a role but we consider 
only one field at a time. 

For conformally invariant models, the scaled gaps Xh(0, L) tend to the 
magnetic scaling dimension for L — > oo ||14||. For nonzero but small t, we 



expect finite-size scaling behavior according to 

X h (t, L)=X h + l/(27r)tL^K 1 (t, l)/dt] t=0 + ■■■ (3) 

where y t the renormalization exponent associated with t. If the temperature- 
like field is irrelevant (y t < 0), X h (t, L) converges to X h with increasing L. If 
it is relevant, Xh(t, L) diverges but for small t Eq. (|]) may still give a good 
description of the finite-size data within the range of available L values. This 
formula serves well for the analysis of the finite-size results Xh(t, L) in terms 
of the phase diagram. 



3 Transfer matrix analysis 

We restrict our parameter space by setting x = y. A preliminary scan of the 
phase diagram in the resulting plane, based on the leading eigenvalues and 
eigenvectors for small L, confirmed the existence of three distinct phases: 

1. A high-temperature or vacuum phase for small x and z. It corresponds 
with the paramagnetic phase in the O(n) spin language. The vacuum 
state is separated by a gap from non-empty states. 

2. A moderately dense phase for large x but not too large z. Type-?/ ver- 



tices are abundant, suggesting Ising disorder as described in Ref. [10 
This is the gapless low-temperature phase of the O(n) spin model. 
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3. At sufficiently large z we observe an even denser phase, dominated 
by type-z vertices. Type-?/ vertices are scarce, suggesting Ising order. 
This phase is of the low-temperature O(n) type, but there exists a gap 
between the leading eigenvalues of the even and odd transfer matrices 
that does not vanish when L — > oo. This is because, on lattices with 
even L, loop configurations composed of only type-z vertices, occur 
only in the case of the even transfer matrix. The odd transfer matrix 
is unable to cover an even L lattice completely with type-z vertices. 

In order to locate lines of phase transitions separating these three regions, we 
set z at some fixed value, and adjust x to its transition value x c . We rewrite 
Eq. (H|) as a function of x, where we may e.g. take t = x — x c in first order to 
represent the relevant temperature field that vanishes at the transition. This 
leads to 

X h (x,L) =X h (x,L + 2) (4) 

where we neglect corrections due to irrelevant fields. Thus numerical solution 
of x of Eq. (U) yields an estimate of x c . Solutions were obtained for even 
system sizes L = 2 • • • 10 in Eq. (f|). We ignored the odd systems because 
they display additional effects due to odd-even alternation. 

Figures [TJ and illustrate how these x c estimates are obtained. They show 
the scaled gaps for various system sizes as a function of x, for z = 0.3 (Fig. 
|T]) and for z = 0.8 (Fig. |2|). Intersections of the curves for successive values 
of L yield estimates of x c . The intersections do not coincide as a consequence 
of corrections to scaling, but seem to converge well with increasing size L. 
The expected corrections behave as negative powers of L and extrapolation 
of the x c estimates was done accordingly. The results are shown in Fig. |3|. 

Another way to investigate the phase diagram is to fix x = y at some 
value and to vary z. For small values of x the leading eigenvalues of the even 
transfer matrix display intersections, leading to a kink in the free energy when 
z is increased: the first-order phase transition between the vacuum phase and 
the dense phase, as already expected in Ref . ]TT| . The intersections provide a 



convenient way to solve for the location of the first-order line: one just solves 
for the value of z where the two largest eigenvalues become degenerate. This 
procedure leads to results that seem to converge well; they are included in 
Fig. H The data shown in this figure are in a full (but, of course, only 
qualitative) agreement with the conjectured phase diagram of Ref. [JR] . 
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The solutions of Eq. (^) do not only yield estimates of the critical point, 
but also for the corresponding magnetic scaling dimension Xh which relates 
to the exponent rj describing the decay of the magnetic correlation function 
as rj = 2Xh- This allows us to verify the expected universal behavior along 
the critical lines. Indeed, as is already illustrated in Fig. [I], the finite-size 
estimates of Xh for the critical line at z < 0.4 are in a good agreement with 
the exactly known value Xh = 5/48 = 0.104166 ■ ■ • for the ordinary O(0) 
or SAW critical point. Furthermore, for z > 0.6 the Xh estimates seem to 



converge well to the expected value |10| Xh = —1/16 for branch 4, which 



is the sum of the low-temperature magnetic exponent (-3/16) and the Ising 
exponent (1/8). The crossover between these different ranges of z appears 
to take place near z = 0.54. At this value, the Xh estimates behave in a way 
consistent with Xh = —5/48 as found for the n = point of branch 3 [|i~l|j . 

These results confirm that the unusual universal behavior reported for 
branches 3 and 4 |L0| is not restricted to the special points that admit exact 
solutions, but apply to larger parameter spaces, at least for n = 0. In the 
x = y plane we locate one point corresponding with branch 3, and a line of 
points that display branch-4 like behavior. Since the n = points of branches 
3 and 4 do not lie in the x = y plane, it appears that the universality class 
of branch 3 covers a 1-dimensional manifold of the x, y, z parameter space, 
and that the universality class of branch 4 covers a 2-dimensional manifold. 
It seems reasonable to expect that a similar picture exists for n ^ 0, at least 
for a range of n values. This will be the subject of further investigations. 

The calculations reported here are modest in size: they required only a 
few hours of CPU time on a Silicon Graphics workstation and a few hours 
on a personal computer. 
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Figure 1: Scaled magnetic g function of the bond weight x at constant 

z = 0.3. Different curves apply to different finite sizes L — 2, 4, 6, 8, 10, and 
12. The steeper the line, the larger L. The scaled gaps at the intersections 
are seen to approach the expected value X h = 5/48 (dotted line). 
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Figure 2: Scaled magnetic g function of the bond weight x at constant 

z = 0.8. Different curves apply to different finite sizes L — 2, 4, 6, 8, 10, and 
12. The steeper the line, the larger L. The scaled gaps at the intersections 
are seen to approach the expected value X h = —1/16 (dotted line). 
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Figure 3: Phase diagram of the 0(0) model for x = y. The vertical scale 
parametrizes the strength of the attractions between loop segments. The 
horizontal scale indicates the weight of vertices visited once by a loop. The 
data points were numerically obtained; their errors are much smaller than 
the symbol sizes. The curves serve only as a guide to the eye. 
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